Doping dependence of antiferromagnetism in models of the pnictides 
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We study the doping dependence of the spin-density- wave (SDW) state in four models of the 1111 
pnictides. The random-phase approximation is used to determine the ordering temperature and the 
ordering vector as functions of doping, and to evaluate the contribution of the various orbitals to the 
SDW instability. In addition to the usual assumption of orbitally rotation-invariant interactions, 
we consider the effect of reduced interactions involving the xy orbital, which are anticipated by 
crystal-structure considerations. We find that changing the relative strength of the interaction in 
the xy orbital tunes the system between different nesting instabilities leading to similar SDW order. 
We identify two models as showing reasonable agreement with experiments, while the other two 
display significant discrepancies, and discuss the underlying differences between the models. 

PACS numbers: 74.50.-|-r, 74.20.Rp 



I. INTRODUCTION 

The discovery of superconductivity in the iron pnic- 
tides has triggered intense research activity striving to 
understand this diverse class of materials. Two pnictide 
families have attracted particular interest due to their 
high critical temperatures, namely the 1111 compounds^ 
i?FeAsO and the 122 compounds^ AFe2As2 {R and A 
are rare-earth and alkaline-earth elements, respectively). 
Both system are antiferromagnets at parent-compound 
filling, and must be chemically doped or subjected to 
significant pressure to become superconducting. The 
close proximity of antiferromagnetic (AFM) and super- 
conducting states in the phase diagram evidences an in- 
timate connection between the magnetic and the super- 
conducting phases.— 1^ 

There is considerable experimental evidence that the 
1111 and the 122 families display a metallic spin-density- 
wave (SDW) state. ^"-'^^ Furthermore, spectroscopic mea- 
surements indicate at most intermediate correlation 
strengths. ^^'^■^ This suggests that the microscopic origin 
of the AFM phase can be understood within a weak- 
coupling picture, which is consistent with ah initio calcu- 
lations identifying the nesting between the electron and 
hole Fermi pockets as responsible for the SDW>i^""— 

The nesting picture of the magnetism implies a strong 
doping dependence of the AFM phase. It is well es- 
tablished that the SDW critical temperature Tc of the 
1111 compounds decreases upon electron doping.^ In 
contrast, very few experiments with hole-doped 1111 
compounds have been reported, and the magnetic be- 
havior remains obscureJ^^— The resistivity anomaly 
characteristic of the onset of SDW in the parent 
compounds becomes much less pronounced upon hole 
doping in Lai„2;Sr:rFeAsO,l'^'^® Pri^^rSr^rFeAsO,!^ and 
Pri_a;Ca2^FeAsO,^*' but it is unclear whether this corre- 
sponds to a suppression of magnetic order. A more ac- 
curate probe of the antiferromagnetism is given by /iSR 
measurements, which in Pri-j^Sr^^FeAsO suggests that 
mesoscopic phase separation allows a substantial frac- 
tion of the system to remain antiferromagnetic up to 



X = 0.2, although it is not known if this is intrinsic— 
The experimental situation is much clearer in the 122 
family, e.g., e/ectron-doping rapidly suppresses the SDW 
in Ba(Fe;^_^Coa;)2As2, the SDW is already absent for 
X = 0.05r^ while for the hole-doped Bai_2:Ka;Fe2As2, 
and the SDW Tc decreases more slowly and vanishes 
at 0.3 < X < 0.4, corresponding to 0.15-0.2 holes per 
single-Fe unit celli^i^ It is interesting to see whether 
this asymmetric doping dependence of the AFM order 
can be replicated by theoretical models. 

The complicated structure of the Fermi surface of the 
pnictides, which consists of two electron pockets and up 
to three hole pockets of various orbital character, allows 
for a number of different nesting instabilities that may 
compete with or reinforce each other. Previous works 
have identified two instabilities that may be responsi- 
ble for the characteristic (vr, 0) order: Many authors^^"— 
identify the most important nesting as that between the 
hole pockets around the F point and the electron pocket 
at the X point of the single-Fe Brillouin zone, which im- 
plies a dominant role of the Fe Sd xz and yz orbitals 
in the formation of the SDW state. On the other hand, 
the role of nesting between the pockets at the Y and 
M points has also been emphasized r^"— where the key 
contribution comes from the Fe 3ci xy orbital. The two 
nesting instabilities are not mutually exclusive, and ap- 
pear at similar doping levels. However, the number and 
location of the hole pockets involved in the magnetic in- 
stability nevertheless has important consequences for the 
spin wave spectrum^i^ and the realized commensurate 
SDW order4i One of the goals of this work is therefore 
to examine the relative importance of the two nesting 
instabilties across the phase diagram. 

The most popular starting point for the theoretical de- 
scription of the pnictides are extended Hubbard models 
possessing orbital degrees of freedom. These models typi- 
cally include up to five of the Fe 3d orbitals^e^iSS^^-i? 
or also the As Ap orbitals42ii^ Since the former dom- 
inate the Fermi surface , ^"^'^"i^^ models constructed in 
terms of the five Fe 3d orbitals are expected to give a 
good account of the physics while keeping the parame- 
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ter space manageable. In the study of these models it is 
common to assume orbitally rotation-invariant interac- 
tions.^^^-^>^-i4.^^^^ This symmetry is however 
broken by the underlying crystal structure of the FeAs 
lattice. Furthermore, upon integrating out the As Ap or- 
bitals to obtain an effective Fe 3(i model, the maximally 
localized Wannier functions of the 3c? orbitals acquire a 
significant As 4p component, and are consequently ex- 
panded compared to the atomic limit The xy or- 
bital has the greatest overlap with the relevant 4p states 
and hence the greatest extent, which implies a weaker 
intraorbital interaction strength. Conversely, the x'^ — 
and — orbitals have only small overlap with the 4p 
orbitals and are thus expected to remain rather compact, 
with stronger interaction potentials. The degenerate xz 
and yz orbitals lie in between these extremes. Miyake et 
al.— have confirmed the expected hierarchy of interaction 
strengths within a constrained random-phase approxima- 
tion (cRPA). The xz, yz, and xy orbitals make up the 
bulk of the Fermi surface?^ and so they are most impor- 
tant when considering the origin of the SDW. The rela- 
tively weaker interaction strength on the xy orbital might 
therefore be expected to play a significant role in selecting 
the leading nesting instability. Investigating this aspect 
of the physics is the second major goal of our paper. 

In order to answer these questions we consider the dop- 
ing dependence of the magnetic order in four different 
five-orbital models of LaFeAsO, proposed by Kuroki et 
al-r^ Graser et al.i^? Ikeda et al.^ and Calderon et al^ 
All four models are based on ah initio band structures. 
Kuroki et al. and Ikeda et al. obtain tight-binding mod- 
els from LDA calculations employing maximally local- 
ized Wannier functions, while Graser et al. fit a Slater- 
Koster tight-binding model to a GGA band structure. 
Calderon et al. propose a Slater-Koster model contain- 
ing a limited number of free parameters, which were 
chosen to best reproduce the LDA band structure. We 
employ the random-phase approximation (RPA) to cal- 
culate the static spin susceptibility in the paramagnetic 
phase, and examine it for divergences as the tempera- 
ture and doping are varied. This allows us to determine 
the limits of the paramagnetic state in an unbiased way 
as we are able to identify ordered states with large unit 
cells that would not be accessible by the usual mean-field 
approaches i^^i'^^'^'^i'^'^'^^i^^ Further insight into the nesting 
mechanisms in the models is obtained by decomposing 
the paramagnetic RPA spin susceptibility into its orbital 
components close to the critical temperature. The com- 
parison of the magnetic phase diagrams of these models 
is the final major goal of our work. Not only does this 
allow us to better understand the mechanisms for antifer- 
romagnetism in the pnictides, but also it helps to identify 
the most realistic models for this system. 

The paper is organized as follows. In Sec.|ll]we intro- 
duce the model Hamiltonian and outline the calculation 
of the RPA spin susceptibility. We proceed to conduct 
a systematic analysis of the four different models in Sec. 
mil This is followed in Sec. IIVI bv a discussion of our re- 



sults and the implications for the understanding of the 
pnictides. We conclude with a summary in Sec. El 



II. THEORY 
A. Model 

We start with a tight binding Hamiltonian that de- 
scribes the non-interacting five-orbital system 



(1) 
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where dl^ ^ ^ {du,u.a) is the creation (annihilation) opera- 
tor for a spin a electron of momentum k in the orbital v. 
The values of for the different models are provided 
in Refs. IH, [H, [s^, and Keeping only local terms in 
the interaction Hamiltonian, we have 

i I'f^ aa' 



"i,i/.(T"i.i/,o-'"i:A':0'"i,M:0' 



(2) 



where the index i stands for the lattice site, v and fi 
stand for the orbitals, and a and a' denote the spin. 
Equation © is usually^e^^M-^dz-ii^^^ simplified 
by assuming orbital-independent interactions for which 
we then have 

Hint = ^ X! X! '^l:''.t"'i,i',i 
i I' 

+ i(2[/ - 5J) ^ ^ ^ ni^,,,ni,^,,, 

i l^^P^ (7,(7' 

i v^p, 

+ '^EE<-t<.i^i.MCfi.Mt- (3) 
i i^T^M 

This is obtained from Eq. ([2]) by setting J/^j^ = U, 
J^f, = J and U„f, = j{2U - 5 J) if ^ fi. The latter 
choice implies invariance of the interaction Hamiltonian 
under rotations in orbital space»^i^ The spin operators 
are expressed in terms of the creation and annihilation 

operators as S^^,, = 5 Z^o-o-' 4,iy,<T<''o-o-"^i,i'o-' , where a is 
the vector of Pauli matrices. The ratio between the local 
Coulomb interaction and the Hund's rule coupling is set 
to U / J = 4. For our calculations we adopt the standard 
assumption that doping does not change the interaction 
and band structure parameters in the Hamiltonian. 

Although the invariance of the interaction Hamil- 
tonian under orbital rotations is a common assump- 
tion^>^^-^i2-ii>H^^>S in general we expect differ- 
ent interaction strengths U^^^, and J^.j/ for inequivalent 
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choices of the orbitals /i, i^. In particular, of the three or- 
bitals which dominate the electronic structure near the 
Fermi surface, ab initio calculations of Miyake et al^ 
predict a much weaker interaction on the xy orbital com- 
pared to the equivalent xz and yz orbitals. In order to 
capture this aspect of the physics we renormalize every 
interaction strength in Eq. ([3]) involving the xy orbital by 
a multiplicative factor of Vxy < 1; all other interactions 
are left unchanged. 



B. Method 

The static paramagnetic spin susceptibility contains 
all necessary information about the magnetic instabili- 
ties of the paramagnetic state. In particular, it diverges 
at the ordering vector as one approaches the critical tem- 
perature Tc- Hence, a temperature-vs. -doping phase di- 
agram for the boundaries of the paramagnetic phase can 
be obtained by determining the highest temperature for 
which l/xs(q, 0) vanishes as the filling is varied, where 
the corresponding q is the ordering vector Q. We obtain 
the static susceptibility from the analytic continuation 
iuJn u> + iO+ of the total susceptibility in Matsubara- 
frequency space and subsequently taking the limit w — > 0. 
The total spin susceptibility is defined as 



Xs(q, iw, 
1 



E 



E 



dre'""^(r.5^(q,T)5^(-q,0) 



(4) 



The Fourier-transformed spin operator 5'^(q) is related 
to the spin operator at site i by 



Due to rotational symmetry in the paramagnetic phase, 
the spin susceptibility can be expressed in terms of its 
transverse part 



X.(q,0) = -x-+(q,0). 
The transverse susceptibility is written as 
X~+(q,iwn) 



I/, /J- 



(5) 



(6) 



where we introduce the generalized susceptibility 
1 ^ rP 



N ^ 

k,k' 



(T^-4+q,..iW'^k,.',t(r)4'-q,p,t(0)^k',^',;(0)).(7) 



The orbitally resolved susceptibilities x 



of particular interest, as they contain information about 
the contribution of the different orbitals to the instability. 
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FIG. 1. (color online) Band structure along hi gh-s ymmetry 



We calculate using the RPA. Summing up the 

ladder diagrams, we obtain the Dyson equation 



where the non-zero elements of Vap-yS are given by 



(8) 



- (1 - V„j)Sa,.,.y] (9a) 

- {I - V„j){Sa,-xy + Sb.,,y)] (9b) 

2J)[1 - (1 - V,y){Sa.,o.y + Sb..y)] (9c) 

- {I ~ Va,y){Sa,-xy + Sb.,,y)] (M) 

with a ^ b. The bare susceptibilities Xiy^'fl' are given by 



Vaaaa = U[l 
Vaabb = J[l 
Vabba - (U - 
Vabab = J[l 



- + (0) / • N 



1 



u,^^'{k)ul^{k)u,>^^,> (k + q)ii*/,^(k + q) 



k s.s' 

nF{Es,k) - 7lF(-Bs',k+q) 



E, 



s,k 



E 



s',k+q " 



(10) 



where np{E) is the Fermi function, E^ are the eigen- 
values of Hq, and Us^^{]i) are the coefficients that trans- 
form the annihilation operators of the diagonalizing ba- 
sis 7s,k into the orbital basis, i.e., d^,]^ = 'u-s,i^i'k.)js,'k- 
We adopt the common approximation to ignore Hartree 
shifts r^i^i^i^ as we assume them to be included in the 
ab initio calculations. We have verified that the inclusion 
of Hartree shifts only leads to small quantitative changes 
in our results. 



III. RESULTS 

In this section we present our analysis of the four dif- 
ferent five-orbital models.— i^i^i^ We first investigate 
the temperature-vs. -doping phase diagrams, calculated 
as described in section IIIBI The summation over k in 
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FIG. 2. (color online) Density of states at the Fermi level as a 
function of the doping S, related by the electron concentration 
hy n — 6 + 5, for the four band structures from Fig. [1] 




FIG. 3. (color online) Critical temperature as a function of 
doping 5 for the four models from Fig. [T] with Vxy = 1, see 
text. Solid lines denote (tt, 0) order, whereas dashed lines 
denote incommensurate order. The interaction strengths were 
chosen as U = 0.875 eV for Kuroki et al.^ U = 0.885 eV 
for Ikeda et al.,^^ U = 1.2232 eV for Graser et al^ and 
U = 1.383 eV for Calderon et al^ 



Eq. pU)) was performed using a 400 x 400 k-point mesh. 
We have checked that a 200 x 200 k-point mesh results 
in only small quantitative differences and so finite size 
effects are expected to be negligible. We ignore ordered 
states with critical temperatures less than lOK. Through- 
out we measure the doping 6 relative to parent compound 
filling, i.e., the electron concentration is n = 6 -I- ^. 

We start by discussing the non-interacting system. In 
Fig. [1] we plot the non-interacting band structure of the 
undoped models close to the Fermi energy. The models 
of Kuroki et o/4^ and Ikeda et al^ are hardly distin- 
guishable within about 0.1 eV of the Fermi surface, and 



so we expect a very similar phase diagram for these two 
models at weak doping. Near to the F point, the Fermi 
surfaces for the model of Graser et al^ almost coincide 
with those of Kuroki et al. and Ikeda et ai, but else- 
where there are significant differences: The hole pocket 
at the M point is very small compared to these models, 
the electron pocket at the X point is also smaller, and 
the {3z^ — r^)-derived flat band at the M point lies much 
closer to the Fermi energy. The model of Calderon et al^ 
is quite distinct from the other models, with highly ellip- 
tical electron pockets, almost degenerate hole pockets at 
the F point, and no flattening of the (Sz^ — r^)-dcrived 
band at the M point. 

The density of states as a function of doping, shown in 
Fig. [21 provides additional insight. The models of Ikeda 
et al.^S. and Kuroki et al^ both show a peak in the den- 
sity of states close io8 = —0.2, although it is significantly 
larger in the latter. This peak arises from the fiat bottom 
of the electron band at the X point. The model of Graser 
et al^ also shows a very high density of states below 
5 = —0.3 that is connected to the (3z^ — r^)-derived flat 
band at the M point. These peaks raise the possibility of 
a competition between AFM and ferromagnetic order at 
strong hole doping. In contrast, the model of Calderon 
et ah'^- has an almost featureless density of states, due 
to the much lower bottom of the electron bands and the 
absence of the (Bz^— r^)-derived flat band at the M point. 

In order to study the doping dependence of the AFM 
phase we first choose the interaction strength U for each 
model such that at 5 = the critical temperature is 
Tc ~ 150 K, close to the experimentally observed order- 
ing temperature. The resulting temperature-vs. -doping 
phase diagrams for all four models are shown in Fig. [31 
All models show an enhancement of Tc for non-zero dop- 
ing. The dome structure of with a maximum at mod- 
erate hole doping [5 ~ —0.1) for the models of Kuroki 
et al^'^ and Ikeda et al^ is consistent the phase dia- 
gram determined by Ikeda et al. in Ref. [s^ . In contrast, 
the monotonous increase of Tc in the model of Graser et 
al^ for strong hole doping, and the enhancement of mag- 
netic order for electron doping in the model of Calderon 
et al.^ are contradicted by experiments. The model of 
Calderon et al. further deviates from experimental find- 
ings by displaying a highly incommensurate ordering vec- 
tor Q = (tt, 0.247r) at zero doping. 

The models of Kuroki et al. and Ikeda et al. are 
most consistent with the reported asymmetric doping- 
dependence of the AFM in the pnictides, although the 
optimal critical temperature is much too high. Note 
that there is no reason to assume that the undoped 
f/iree-dimensional parent compounds are best described 
by choosing precisely the value (5 = in these two- 
dimensional models. Moreover, we will see that a mod- 
erate increase of the onsite energy of the xy orbital shifts 
the peak position to (5 = for some of the models. In- 
stead of specifying Tc at (5 = as in Fig. [31 it is therefore 
more reasonable to search for an interaction strength giv- 
ing a dome shape of Tc and an ordering temperature of 
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FIG. 4. (color online) (a) Doping dependence of the critical 
temperature for the model of Kuroki et alM- at different values 
of V^y. To maintain r°P' « 165 K, we choose U = 0.8 eV for 



1, f/ = 0.845 eV for K,„ = 0.9, U = 0.885 eV for K,. 



0.8, and U = 0.915 eV for 14^ = 0.7. (b) Doping dependence 
of the critical temperature for different shifts of the xy-orhital 
on-site energy Ssxy with respect to the value given in Ref. 42. 
For constant optimal doping critical temperature T°'^^ « 165 
K we choose U = 0.8 eV for Ss^y = eV, C/ = 0.782 eV 
for Se^y = 1.04 eV, U = 0.78 eV for Ss^y = 1.09 eV, and 
U = 0.785 eV for Se^y = 1.13 eV. 



jiopt ^ K at a non-zero optimal doping S. The lower 
jiopt sjjQ^ifj g^igQ allow us to unambiguously identify the 
leading nesting instabilities responsible for the SDW. Be- 
low we construct the phase diagram of each model ac- 
cording to this argument. 



A. The model of Kuroki et al. 

The phase diagram for the model of Kuroki et al^ is 
shown in Fig. 21 Initially focusing on the case of Vxy = 1, 
we still find a dome of (tt, 0) order centered at the optimal 
doping 5 — —0.093, but there is no magnetic order at 
(5 = 0. In addition to the AFM dome, there is also a 
region of ferromagnetic order around 5 = —0.2. The 
highest critical temperature of the ferromagnetic state 
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FIG. 5. (color online) Fermi surface for the model of Kuroki 
et alM' (heavy lines) superimposed with the Fermi surface 
shifted by (a) (vr.O) for 5 = -0.093 and (b) (0.957r, 0) for 
5 — —0.105 (thin lines), corresponding to the ordering at 
optimal doping for Vxy ~ 1 and Vxy = 0.7, respectively. 



1000 



-r 100 



OS 




(0.93t,0) (11,0) (71,0. l7t) 
a) 1 



(b) 



(0.97t,0) (71,0) (71,0. 17I) 



1000 



^ 100 r 




10 " .-J^ \«.\ 



1000 



100 



10 







,22 
yz,3z -r 






xz,xz ■ ^— 
yz,yz — — 
X y , X y - - ■ 

^ yz,xz — ■ 






total 


i*^ * 
'/ 





(0.971,0) (7t,0) (7t,0.l7t) 

c) q 



(0.971,0) (71,0) (Jt,0.1Jt) 



(d) 



FIG. 6. (color online) Total susceptibilities and orbitally re- 
solved contributions Xutj" at T = 180 K for the model of 
Kuroki et aZ.— at (a) Vxy = I, S = -0.093, (b) Vxy = 0.9, 
S = -0.092, (c) Vxy = 0.8, S = -0.096, and (d) Vxy = 0.7, 
S = -0.105. 



occurs close to the peak in the density of states, see Fig. 
[21 consistent with the Stoner criterion. 

The mechanism responsible for the (tt, 0) order can be 
observed most clearly at the optimal doping 6 = —0.093. 
As shown in Fig. [Sja), here we find excellent nesting of 
the a;?/-orbital-dominated parts of the electron pocket at 
the Y point with the xy-derived hole pocket at the M 
point, suggesting that these Fermi surfaces play the lead- 
ing role in the AFM instability. A further hint that the xy 
orbital is most important for the AFM order comes from 
the observation that the optimal doping can be shifted to 
i5 = by increasing the on-site energy of the xy orbital 
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FIG. 7. (color online) Doping dependence of the critical 
temperature for the model of Kuroki et alM- for Vxy ~ 0.7 
(dashed) and for the fully orbital-dependent interactions pre- 
dicted in Ref. [Z^ (solid). For the latter, the interactions 
strengths were scaled by 0.406 to find a maximal ordering 
temperature of Tc ~ 165 K at <5 > —0.2. Both calculations 
yield an ordering vector Q — (0.957r, 0) at the maximum Tc 
close to 5 = —0.1. 



face. By changing Vxy we can therefore select the domi- 
nant nesting instability of the system. 

Our treatment of the weaker interaction on the xy or- 
bital neglects the likely different interaction strengths 
involving the other inequivalent orbitals. To test our 
approximation, we compare the Vxy — 0.7 phase dia- 
gram of Fig. m^a) with the phase diagram calculated us- 
ing Miyake et al.'f^- cRPA values for U^fj_ and J^^ in 
Eq. As shown in Fig. [71 upon suitably rescaling 

the cRPA interaction potentials, we find excellent agree- 
ment between the two phase diagrams for the physically 
reasonable doping regime 5 > —0.2. We note that in 
the fully orbital-dependent results we have to choose a 
slightly larger value of Uyz,yz than in our Vxy = 0.7 
calculations in order to achieve r°P* ~ 165 K for the 
peak around 5 ~ —0.1; the origin of this discrepancy is 
likely the weaker exchange interaction J^^ < 0.25 U^i, 
predicted by Ref. \^ This also indicates that the AFM 
order does not crucially depend upon the ratio U / J in 
the weak-coupling regime. For strong doping, larger de- 
viations appear, in particular there is no ferromagnetism 
for the cRPA interactions, and the critical temperature 
of the incommensurate AFM state at (5 ~ —0.4 is higher 
although the ordering vector is similar. 



by ^ 0.1 eV [see Fig. HJ^b)], and always coincides with 
good nesting of the Y and M pockets. 

To examine the effect of the likely weaker interactions 
in the xy orbital, in Fig. |4l[a) we plot the evolution 
of the phase diagram upon reducing Vxy while simulta- 
neously increasing U such that T°p' remains constant. 
At Vxy =0.7, which is close to the value predicted by 
Miyake et al.^ we find that the optimal doping of the 
AFM dome shifts to 5 = —0.105, and the ordering vec- 
tor at this filling becomes weakly incommensurate with 
Q = (0.957r,0). The Fermi surface plotted in Fig. E^b) 
shows that this corresponds to good nesting between the 
hole pockets at the F point and the electron pocket at 
the X point, where the best-nested segments of the Fermi 
surface have mostly yz-orbital character. The necessary 
increase of /7 by ~ 10% when Vxy is reduced implies 
that the AFM instability due to F-X nesting requires a 
significantly higher interaction than the Y-M nesting to 
produce a realistic r°P*. Additional magnetically ordered 
states appear at strong doping 5 « —0.42 with optimal 
ordering vector Q = (0.447r,0). At these doping levels, 
however, our assumption of rigid bands is questionable 
and so the physical relevance of these results is doubtful. 

Greater insight into the origin of the AFM order can be 
achieved by examining the orbitally resolved susceptibil- 
ities at optimal doping just above Tc, see Fig. [51 As Vxy 
is decreased, the dominant contribution shifts from the 
susceptibilities xZi^ involving the xy orbital to those in- 
volving the yz orbital, and the peak in moves from 
(tt, 0) to an incommensurate vector. This is in perfect 
agreement with the observed change of the nesting from 
xy-dominated to j/z-dominated parts of the Fermi sur- 



B. The model of Graser et al. 

The phase diagram for the model of Graser et al^ 
is shown in Fig. EJa). At Vxy = 1 we find a dome of 
commensurate AFM order centered at i5 = —0.12 and 
also a small incommensurate dome with Q = (0.957r, 0) 
aX 5 K, —0.05. At 5 < —0.2 we observe an incommen- 
surate AFM state with ordering vector Q = (0.7l7r, 0), 
and a high critical temperature, which strongly increases 
with hole doping. The strong tendency to AFM order at 
strong hole doping occurs only in this model, and is likely 
connected with the (Sz^ — r^)-derived flat band at the M 
point. Indeed, a Fermi surface due to this band appears 
at the critical doping level for the incommensurate AFM 
order. Furthermore, the ordering at strong doping can 
be suppressed by increasing the on-site energy of the xy 
orbital, which effectively lowers the flat band at the M 
point. This also shifts the optimal doping of the (tt, 0) 
dome at (5 = —0.12 towards zero, again suggesting an 
important role for the xy orbital in the (tt, 0) order. 

Focusing our attention on the regime of moderate dop- 
ing, —0.2 < (5 < 0, we find that as Vxy is reduced, 
the small incommensurate dome grows and becomes the 
leading instability at Vxy = 0.7. At this value of Vxy, 
the ordering vector at optimal doping 5 = —0.06 is 
Q = (0.937r,0). The orbitally resolved susceptibilities 
at the optimal doping of the (tt, 0) order at Vxy — 1 and 
the incommensurate order for Vxy — 0.7 are plotted in 
FigS-EUb) and (c), respectively, while the nesting of the 
corresponding Fermi surfaces are shown in Figs.[5Jd) and 
(e). These results are very similar to those obtained for 
the model of Kuroki et al.^ see Figs. [S] and [51 We hence 
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FIG. 8. (color online) (a) Doping dependence of the critical 
temperature for the model of Graser et al^ at different values 
of V^y. To maintain r°P* ^ 165 K, we choose U = 1.09 eV 
for = 1, ?7 = 1.146 eV for V^^y = 0.9, U = 1.191 eV for 
V^y = 0.8, and [/ = 1.222 eV for = 0.7. (b) The total 
susceptibility and the largest Xui^ ^or Vxy = 1, T — 180 K, 
and the optimal doping 5 = —0.137 of the (7r,0) state, (c) 
Same as in (b) but for Vxy = 0.7, T = 180 K, and the optimal 
doping 5 = —0.065 of the incommensurate state, (d) Fermi 
surface at (5 = —0.137 (heavy lines) with the same Fermi 
surface shifted by (vr, 0) superimposed (thin lines), (e) Fermi 
surface at 5 = —0.065 (heavy lines) with the same Fermi 
surface shifted by (0.937r, 0) superimposed (thin lines). 



FIG. 9. (color online) (a) Doping dependence of the critical 
temperature for the model of Ikeda et al^ at different values 
of Vxy To maintain r°P* ^ 165 K, we choose U = 0.808 eV 
for Vxy = 1,U = 0.852 eV for Vxy = 0.9, U = 0.892 eV for 
Vxy = 0.8, and U = 0.925 eV for Vxy = 0.7. (b) The total 
susceptibility and the largest Xv/i for '^xy = 0.7, T — 160 K, 
and S = -0.12. (c) Same as in (b) but at 5 = -0.06. (d) 
Fermi surface a.t 5 — —0.12 (heavy lines) with the same Fermi 
surface shifted by (0.957r, 0) superimposed (thin lines), (e) 
Fermi surface at 5 = —0.06 (heavy lines) with the same Fermi 
surface shifted by (7r,0) superimposed (thin lines). 



conclude that again a reduction of V^^y tunes the system 
from the a;?/-dominated to the yz-dominated instability, 
although the optimal dopings for the two instabilities are 
more widely separated than in the model of Kuroki et al. 
This can be explained by the small M pocket at zero dop- 
ing in Graser et al.'s model,-^'^ which implies a much larger 
change of the doping in order to optimize the Y-M nest- 
ing. Furthermore, due to the smaller size of the electron 



pockets at parent compound filling, the incommensurate 
AFM phase is stabilized at weaker doping. 



C. The model of Ikeda et al. 

The phase diagram for the model of Ikeda et al~^ at 
Vxy = 1 [Fig. m^a)] is very similar to that for the model 
of Kuroki et al. ^ although ferromagnetism is not found 
near 6 — —0.2, consistent with the smaller peak in the 
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density of states [Fig. [2]. The (7r,0) order is again dom- 
inated by the xy orbital, and the nesting of the Y and 
M pockets is primarily responsible for the AFM state, in 
agreement with Ref. l39l . Despite the very similar band 
structure to the model of Kuroki et ai, the phase di- 
agram for Vxy < 1 shows a significant difference: The 
(tt, 0) order does not vanish for Vxy = 0.7, but instead 
moves to the low-doping half of the AFM dome while 
the other half is incommensurate. As for Kuroki et al.'s 
model, however, the yz orbital dominates the magnetism 
over the full doping range. This can be seen in two repre- 
sentative plots of Xu^t Figs.EJb) and (c) at 6 = —0.12 
and S — —0.06 which correspond to incommensurate 
Q = (0.957r,0) and commensurate Q = (7r,0) order- 
ing vectors, respectively. As revealed by Figs. El^d) and 
(e), at both of these doping levels there is excellent nest- 
ing of the X and T pockets. The ordering vector tracks 
the continuous evolution of this nesting vector across the 
dome, from slightly incommensurate to commensurate. 
The close similarity of this model to the one of Kuroki et 
al. shows that very small differences in the band structure 
play a major role in determining the ordering vector. 



D. The model of Calderon et al. 

The phase diagram for the model of Calderon et al.^ 
shown in Fig. fTUf a). is in stark contrast to those for the 
other models. The magnetic order is optimized at strong 
electron doping, 5 « 0.18, and occurs at the highly in- 
commensurate ordering vector Q = (7r,0.397r). Although 
this is inconsistent with experimental findings for the 
pnictide systems, it is nevertheless interesting to examine 
the origin of this AFM state. An important clue comes 
from observing that the phase diagram hardly changes 
when Vxy is decreased, although the interaction U has to 
be increased in order to keep T°p* constant. This indi- 
cates that the xy orbital is almost exclusively responsible 
for the magnetic ordering, which is confirmed by the or- 
bitally resolved susceptibilities presented in Figs. fTOlb) 
and (c). The incommensurate ordering vector at optimal 
doping gives excellent nesting of the xy pocket at the M 
point and the xy-dominated tip of the electron Fermi sur- 
face at the Y point [Fig. [TOld)]. revealing an unexpected 
similarity to the nesting instabilities in the other models. 

The stabilization of AFM order by electron-doping in 
the model of Calderon et al. follows from the observation 
that the magnetic order arises only from Y-M nesting. 
In Fig. [T] it can be seen that at 5 = the circular hole 
pocket at the M point has almost the same radius as 
in the models of Kuroki et al.'^^ and Ikeda et al.^^ In 
contrast, the ellipticity of the electron pockets is highly 
exaggerated in the model of Calderon et al, and indeed 
the minor axis of the electron pockets is much smaller 
than the diameter of the M-point hole pocket. As such, 
it is necessary to raise the chemical potential (i.e., dope 
with electrons) to optimize the nesting between these two 
Fermi surfaces. 
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FIG. 10. (color online) (a) Doping dependence of the critical 
temperature for the model of Calderon et al^ at different 
values of V^^y. To maintain r°P* 165 K, we choose U = 
1.269 eV for V^y = 1, f/ = 1.378 eV for V^y = 0.9, U = 
1.498 eV for V:,y = 0.8 and U = 1.635 eV for = 0.7. (b) 
The total susceptibility and the largest Xvi^ for Vxy = 1, T = 
180 K, and 5 = 0.176. (c) Same as in (b) but for V^y = 0.7. 
(d) Fermi surface at 5 = 0.176 (heavy lines) with the same 
Fermi surface shifted by (tt, 0.397r) superimposed (thin lines). 



IV. DISCUSSION 

Our study of the four different models for the 1111 
pnictides shows a clear distinction between the models 
of Kuroki et al.*^ Graser et al.^ and Ikeda et al^ on 
the one hand and the model of Calderon et al^ on the 
other. In the former, a SDW is stabilized for hole dop- 
ing with commensurate or near-commensurate ordering 
vector, whereas in the latter we find a strongly incom- 
mensurate AFM state upon electron doping. The be- 
havior of Calderon et al.^s model can be explained by 
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the unrealistically high elhpticity of the electron pock- 
ets. Calderon et al^ started from an eight-orbital model 
with nearest-neighbor Fe-Fe and Fe-As hopping and then 
removed the As orbitals within perturbation theory. By 
using the Slater-Koster approach, they were able to write 
the eighteen hopping integrals as functions of only four 
overlap integrals and the Fe-As bond angle; in the other 
models, these eighteen hopping integrals are free parame- 
ters, which gives much greater freedom in fitting the band 
structure. One can speculate that the smaller parameter 
space available to Calderon et al. is responsible for the 
unphysical doping dependence in their model. Among 
the remaining models, the monotonous increase of Tc for 
the model of Graser et al, see Figs. [3] and [HJ a), also con- 
tradicts available experimental data^ii"— i^^i^ The com- 
parably high (3z^ — r^)-derived flat band at the M point is 
responsible for this strong tendency to AFM order.— In 
contrast, the models of Kuroki et al. and Ikeda et al. show 
a dome of (7r,0) order centered at moderate hole dop- 
ing 5 « —0.1; for realistic maximum critical temperature 
Tc ~ 165 K there is no ordering at zero nominal doping, 
(5 = 0. Similar behavior is seen in the model of Graser 
et al. for sufficiently small interaction strength. In the 
first three models there are hence strong (tt, 0) spin fluc- 
tuations at weak doping, consistent with previous stud- 
jgg39j42j43 ^Yie asymmetric doping dependence of the 
SDW phase observed in some experiments i^ii2r^ 

Our results again demonstrate the sensitive depen- 
dence of the magnetic order on small details of the band 
structure in the weak-coupling limit/^ This is illustrated 
by the major discrepancies between the model of Graser 
et al. and the more realistic models of Kuroki et al. and 
Ikeda et al., despite the very similar band structures of 
these three models. It is significant that the models of 
Kuroki et al. and Ikeda et al. were obtained by fitting to 
ah initio results for the experimental crystal structure of 
the FeAs planes, whereas the model of Graser et al. is 
based upon ab initio calculations for a relaxed structure. 
These two structures have rather different Fe-As bond 
angles, which has been identified as a crucial control pa- 
rameter for pnictide physicsi^i^i^ 

In addition to the major role of the band structure, im- 
portant details of the antiferromagnetic order were found 
to be controlled by the interaction strengths involving 
the xy orbital relative to the other orbitals. Under the 
common assumption of orbitally rotation-invariant inter- 
actions, the models of Kuroki et al^ and Ikeda et al^ 
display (tt, 0) order at optimal doping (maximum SDW 
critical temperature Tc). Upon reducing the interactions 
involving the xy orbital, the optimal doping slightly shifts 
and the ordering vector tends to become weakly incom- 
mensurate: At the lowest reduction factor Vxy = 0.7, 
the AFM dome for the model of Kuroki et al. becomes 
entirely incommensurate, while for the model of Ikeda 
et al. there is a continuous change from incommensurate 
to commensurate order as one moves from stronger to 
weaker doping. Reducing Vxy in the model of Graser et 
al.^ we find that the dome with (tt, 0) order is almost 



obscured by a second dome with weakly incommensurate 
order that has a smaller optimal doping. 

We find that in all models commensurate (tt, 0) or- 
der is realized for an extended doping range. This shows 
that commensurate order can be stabilized in purely elec- 
tronic models, but does of course not imply that mech- 
anisms beyond the models considered here are unim- 
portant in this respect. Indeed, the experimentally ob- 
served robustness of commensurate order suggests that 
some additional stabilizing mechanism is required, espe- 
cially in the scenario of reduced Vxy. Accounting for the 
three-dimensionality of the Fermi surface may improve 
the nesting with a commensurate ordering vector. Most 
likely, though, the weak incommensuration we have found 
will result in commensurate order if magneto-elastic cou- 
pling is taken into account. 

Our results allow us to distinguish the dominant nest- 
ing instability in the various models. The (tt, 0) order 
observed for orbitally rotation-invariant interactions orig- 
inates mainly from the good nesting between the xy- 
derived parts of the Y pocket and the M pocket, as 
identified in Refs. IsB - fsol For reduced interactions in- 
volving the xy orbital, however, the nesting instability 
between the F and X pockets gives the highest order- 
ing temperature, which corresponds to the scenario in 
Refs. |25| - [34| . This mechanism is dominated by the yz 
orbital and produces a SDW with the ordering vector 
Q = (Qx, 0), where O.Ott <Qx <t^- 

Our results therefore confirm the two proposals in the 
literature for a nesting instability in the pnictides, with 
their relative importance tuned by the parameter Vxy. 
This implies a crucial role for the As 4p orbitals, as their 
hybridization with the Fe id orbitals is ultimately respon- 
sible for the strong orbital dependence of the interaction 
potentials in an effective 2>d theory.™ It is therefore some- 
what unsatisfying that the As orbitals do not enter the 
calculations more directly. Indeed, keeping the As Ap or- 
bitals reduces the variation in the size of the Fe Zd Wan- 
nier functions, and hence the interaction strength is likely 
to show much less pronounced orbital dependence A 
detailed comparison of the spin fluctuations in a realistic 
3d-4p model with those in a pure 3d model is therefore 
desirable4^ 



V. SUMMARY 

In this work we have presented an analysis of the insta- 
bilities responsible for magnetic order in the 1111 pnic- 
tides. Using the RPA we have determined the doping de- 
pendence of the SDW critical temperature and ordering 
vector in four different five-orbital models. For the three 
models proposed by Kuroki et al.,^ Graser et al.^ and 
Ikeda et al.,^^ we find that the observed (7r,0) magnetic 
order is stabilized for hole doping, while in the model 
of Calderon et al^ an incommensurate SDW phase ap- 
pears at electron doping, contradicting experimental re- 
sults. We have studied the relative importance of the two 
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known nesting instabilities across the phase diagram of 
these models, and have identified the relative interaction 
strength in the xy orbital as a parameter that tunes the 
dominant mechanism leading to magnetic order. We have 
identified two models as giving particularly good agree- 
ment with experiment, and discussed the band structure 
features which lead to the poorer agreement for the oth- 
ers. 
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